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Abstract 

<i ■ 

Q_ | 1 The status of coupling constant unification (with and without a unification 



of Yukawa couplings) is discussed. Uncertainties associated with the input cou- 
<\^\ pling constants, nib an d mt, threshold corrections at the low and high scale, and 

possible nonrenormalizable operators are described and a discrepancy between 
effective and physical scales is pointed out. Theoretical uncertainties in the 
^ ! predictions of a 8 (Mz), rnb, and the unification scale, Mq, are discussed and 

^ estimated. Constraints on the super-partner spectrum are found to be weak if 

uncertainties associated with the high-scale are included. However, requiring 
h b = h T at M G excludes 3 £ tan/3 £ 40 (for m t £ 200 GeV). 

The standard model (SM) couplings were recently shownQ] to meet at a point, 
Mq ~ 10 16 — 10 17 GeV, when extrapolated to high energy assuming a grand desert 
and the spectrum of the minimal supersymmetric standard model (MSSM) above the 
weak scale (a minimal one-step scenario). This is still true if one admits additional 
U\ factors or a small number of complete multiplets of a gauge group of some grand 
unified theory (GUT). (The latter affect the predictions only at the two- loop level.) 
Otherwise, it would be difficult to relax the assumption on either the spectrum or 
the desert scenario without relaxing the other, i.e., introducing additional matter at 
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some intermediate scales. In order to maintain the predictability of the model, we 
assume hereafter a minimal one-step scenario. 

For certain values of the t-quark (pole) mass, m t , and of the two Higgs doublet 
expectation value ratio, tan/5, Yukawa couplings of the third family fermions also 
unify at ~ M G (when extrapolated under the same assumptions) , i.e., hb = h T 
(= h t ), as is implied by certain SU5 (SO 10) and similar GUT's. (One usually assumes 
that some perturbation modifies the couplings or the masses of the two light families 
where, in principle, similar relations should, but do not, hold.) While the coupling 
constant extrapolation is decoupled to a good approximation from that of the Yukawa 
sector, the latter is controlled by the coupling constants. It is the balance between 
the coupling constants and the Yukawa couplings that determines the infra-red fixed 
points in the Yukawa coupling renormalization flow. We will examine the status of 
coupling constant unification first Q. Then, we will further assume Yukawa unification 
at ~ Mq, and use that assumption to constrain the m t — tan/5 plane While mt 
effects can be treated as a correction term in the former case, m t is a free parameter 
in the latter. a s will be fixed by the unification, and as we point out below, is thus a 
quadratic function of m t . 

The naive scenario has to be somewhat relaxed in order to obtain a more realistic 
picture. Below we will perturb the unification and desert assumptions by taking 
m t > Mz, an arbitrary sparticle spectrum (below the TeV scale), and by considering 
a split spectrum and nonrenormalizable operators (NRO's) at the high-scale. For the 
latter, we will assume a minimal SU5 GUT. We do so in order to explicitly realize 
the magnitude of the effects and the role of the constraints coming from proton decay 
at the loop-level 0]; however, this can be easily generalized 0. (Having larger GUT 
gauge groups does not imply larger corrections.) We assume that all the corrections 
are consistent with perturbative treatment, and that there is no conspiracy among 
the different correction terms. We will add theoretical uncertainties in quadrature as 
a guideline only. 

Let us then write 

— -— = — + kt + 6i + Hi - A; for 2 = 1,2, 3, (1) 

ai(M z ) a G 

where we can neglect the two- loop contribution from the Yukawa sector, Hi. All 
other uncertainties and corrections (including conversion to DR) are included in the 
correction functions A$. a G is the coupling at M G and t = ^ In jf- is a convenient 
parametrization of the unification point, Mq. 6, are the respective one-loop /9-function 
coefficients and 9i are the two- loop corrections. One can then get expressions for a@, 
t, and either a s (Mz) or the weak angle, s 2 (Mz), in terms of only a(M z ) and either 
s 2 (Mz) or a s (Mz). Q% are calculated iteratively and Aj determine the theoretical 
uncertainties. 
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s 2 (Mz) is strongly correlated with m t . It is then useful to define a (MS) m t - 
independent quant ity||, Sq = 0.2324 ± 0.0003; i.e., s 2 [M z ) obtained for a fixed m to 
(= 138 GeV), and use Sq to predict a s (M z )- One can then account for m t ^ 138 
GeV (i.e., 113 < m t < 159 GeV from precision electroweak data) by including the 
leading (quadratic) mt dependence in the correction functions Aj||. (This is accurate 
up to small logarithmic corrections). By following the above procedure, (leading) m t 
effects are treated consistently and we bypass the ~ 8% input uncertainty in the 
a s (M z ) = 0.12 ±0.01 range extracted from experiment. The ±0.01 input uncertainty 
would have induced an uncertainty comparable to the theoretical ones if we were 
using a s to predict s 2 . Setting A, = 0, as in the naive calculation, we obtain a^ 1 = 
23.41 ± 0.04, t = 5.30 ± 0.01, and a s (M z ) = 0.125 ± 0.001. Turning on m t > M z 
correlates a s (Mz) with m t , an effect that has to be taken into account when, e.g., 
constraining the m t — tan (3 plane. 

Aside from m t effects, A, consists also of one-loop threshold corrections (which 
is consistent with a two-loop calculation) and NRO effects. The former have to 
be accounted for at both scales @, and analytic treatment is then more instructive 
and convenient. NRO that renormalize (and thus split) the couplings at Mq are 
suppressed by the ratio of Mq ^ 10 16 GeV over some larger scale. Even so, the 
effect can be comparable to two-loop effect s|3|, 0. We parametrize the A$ functions 
in terms of 7 effective parameters (aside from m t ) corresponding to the sparticle and 
Higgs doublet spectrum (Ml, M 2 , M 3 ); the high-scale thresholds (My, M 2 4, M 5 , 
for the heavy components of the vector, adjoint and complex Higgs supermultiplets, 
respectively); and an effective NRO strength, rj (which can have either sign). All 
these are described in detail in Ref. 3. The corrections to the a s (M z ) prediction 
due to each of these parameters are illustrated in Figure 1. The corrections are all 
comparable and have no fixed sign, i.e., 

a s (M z ) « 0.125 ± 0.001 ± 0.005±g;gg| ± 0.002 ± 0.006, (2) 

or a s (M z ) pa 0.125 ± 0.010, which is in good agreement with the data. The uncer- 
tainties quoted in (fj) correspond to those of a and Sg, reasonable choice of ranges 
of sparticle (and Higgs) spectra, high-scale spectra, 113 < m t < 159 GeV, and to 
NRO effects, respectively. They are intended to serve as an order of magnitude esti- 
mate only. (Proton decay constrains M 5 and thus constrains the high-scale threshold 
negative contribution to the uncertainty, unless one turns to simple extentions or 
eliminates the colored triplet Higgs from the spectrum.) 

To further illustrate the underlying formalism, let us discuss in greater detail 
the parameters Mj. They are defined as weighted sums, 

_ bi , M C bf SSM - bf M Mi r . i o q 

> — In — - = — In for i = 1, 2, 3, (3) 

M z 2tt M z ' ' ' V ; 
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where the summation is over all the sparticles and the heavy Higgs doublet and b\ is 
the ^-particle contribution to the respective /3-function. Mj determines the low-scale 
threshold contribution to Aj and it is a straight-forward exercise to calculate it for 
any spectrum^]. (Similar model independent parameters can be defined at the high- 
scale[[J.) The parameters are sensitive to the split between, for example, colored and 
uncolored sparticles. A as is proportional to the combination 25 In jj± — 100 In ||| + 

561n||| (= —19 In As ^ Y ), and is potentially large and negative; e.g., for a spectrum 
degenerate at M S usy- However, in supergravity-inspired MSSM M 3 ~ M gaugino ; 
Mi ~ m , //; and M 2 ~ \i. Thus, in general M 2 ^ M 1 and M 2 ^ M 3 , and the 
correction can be positive, depending on the split. Rather than evaluating A as in 
terms of a common mass, Msusy, we can invert the logic and use the above expression 
to define an effective scale, Asusy, which is the relevant scale in the problem. Asjjsy 
can be as low as a few GeV and does not contain by itself any physical information 
(nor should we expect it to). One should therefore calculate Mi, M 2 and M3, which 
determine also the low-scale corrections to «g and t, and contribute to the corrections 
to mj. It is always possible to define an effective scale in the context of a specific 
prediction, but unlike the Mj parameters, such scales do not have an obvious physical 
interpretation. 

Similarly to (@) we have (using the same ranges for the parameters) 

t » 5.30 ± 0.01 ± 0.09ig;g ± 0.02 ± 0.025, (4) 

which corresponds to 10 16 ^ Mq ^2x 10 17 GeV. t is insensitive to 77 (i.e., to NRO's) 
but increases significantly if the heavy color octet and S'6 r 2-triplet Higgs coming from 
the adjoint representation (M 24: in the above notation) are somewhat lighter than the 
other heavy thresholds. One can then increase the predicted Mq up to ~ 5 x 10 17 
GeV while maintaining a successful prediction of a s by postulating a large (but still 
consistent with a perturbative treatment) and negative rj. (The above considerations 
relax proton-decay constraints 0.) However, it is not straight-forward to realize such 
an interplay in string-inspired GUT-like models that, in general, do not admit adjoint 
and other large representations. 

Finally, we require also h^Mo) = h T (Mo)- There is no expression similar to (|IJ) 
for the (two-loop) Yukawa couplings and one has to turn to numerical integration. 
By solving a set of six two-loop renormalization group equations [[L(J (neglecting the 
two light family Yukawa couplings and flavor mixings) we can predict the 6-quark 
(current) mass, mj, as a function of m t and tan/3. The leading dependences on m t 
are from h t and the a s (M z ) prediction (the s 2 (M z ) input). It is useful to define 
m b = p _1 x mJJ, where m" is calculated to two-loop (numerically) using the ideal 
desert and unification assumptions (aside from m t effects). p _1 is a multiplicative 
correction function (e.g., the equivalent of the additive A as , only that A aa in Eq. 



4 



(0), for example, includes m t corrections and p _1 does not). The correction function, 
p _1 , consists of corrections to a s (M z ), ag, and Mq (t) described previously, as well 
as threshold effects that correct the running of the Yukawa couplings either directly 
(by modifying the respective /3-function) or via the modified running of the coupling 
constants (and in particular, 03). One could also allow a small arbitrary split between 
the two unification points. All these are described in detail in Ref. 4, where we studied 
the effects using (approximate) analytic correction expressions. We obtain pa 
1 ± 0.15, and thus, the constraint 0.85m°(5 GeV) < 4.45 GeV. The 0.85 correction 
factor is estimated for a conservative choice of ranges for the various parameters, and 
is more sensitive to high-scale corrections to the coupling constant unification than 
to the details of the sparticle spectrum (for sparticles below the TeV scale) . We also 
use a conservative estimate of the upper bound on the current mass. 

That constraint, in addition to requiring a perturbative Yukawa sector (h < 3 up 
to Mq), allows only two branches in the m t — tan (3 plane, as illustrated in Figure 
2. The lower and upper bounds on tan/3 and the upper bound on m t are from 
perturbative consistency (i.e., are determined by the fixed-points). The width of 
each branch is determined by the 0.85 correction factor and the current mass upper 
bound. m t correlations increase the upper bound on m t (~ 215 ± 10 GeV), and as a 
result, the allowed area is not sensitive to m t ^ 200 GeV. For m t ^ 200 GeV some 
intermediate values of tan/5 are allowed. Otherwise, tan/? is strongly constrained, 
i.e. 0.6 ^ tan/3 ^ 3 or 40 ^ tan/5 ^ 60. The latter may be further excluded for 
m t ^ 170 ± 10 GeV by requiring radiative breaking of SU2 x U\. 

We can write an explicit constraint on the sparticle (and Higgs doublet) spectrum 
from coupling constant unification, i.e., 



0.110 ^(A, z ) + (a ° (M2))i 



28tt 



Mi , M 2 , M 3 

25 In — i- - 100 In — ± + 56 In 



±0.008 £ 0.130, 



M z M z M z 

(5) 

where a° s here is the prediction for ideal desert and unification (i.e., 0.125). Requiring 
in addition H^Mq) = h T (Ma) implies a similar, but weaker, constraint which depends 
also on max (M gluino , M squark ), min {M gluino , M squark ), and on the heavy Higgs doublet 
mass, Mu. (See Ref. 4.) High-scale effects weaken the constraints significantly. 

To summarize, coupling constant unification agrees very well with the data. Yukawa 
coupling unification strongly constrains tan/? (independent of m t ^ 200 GeV), and 
is in agreement with the data only in a small area of the m t — tan (3 plane. The 
above analysis may be further incorporated, e.g., in a one-loop sparticle spectrum 
analysis [[!]]], and, in principle, can give hints on the structure of the high-scale physics 
(e.g., here we showed that tan/3 pa 3 — 40 will disfavor a large class of models). 
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Figure 1: Contributions of individual correction terms - the SUSY effective mass 
parameters Mf, the heavy thresholds at the high-scale; the t-quark mass; and NRO's 
at the high-scale - to the a s (Mz) prediction. The error bar on the a s (Mz) range ex- 
tracted from experiment (dash-dot line) and the two-loop contribution to the a s (M z ) 
prediction (dotted line) are given for comparison. 



Figure 2: The wF t oe — tan/5 plane is divided into five different regions. Two areas 
(low- and high-tan f3 branches) are consistent with perturbative Yukawa unification 
{h b = h T at Mq) and with 0.85m°(5 GeV) < 4.45 GeV. Between the two branches 
the 6-quark mass is too high. For a too low (high) tan/?, h t (h b ) diverges. The 
strip where all three (third-family) Yukawa couplings unify intersects the allowed 
high-tan (3 branch and is indicated as well (dash-dot line). Corrections to the ht/hb 
ratio induce a ~ ±5% uncertainty in the m? t ole range that corresponds to each of 
the points in the three- Yukawa strip. a s (Mz), «g, and the unification scale used in 
the calculation are the ones predicted by the MSSM coupling constant unification, 
and are sensitive to the t-quark pole mass, mf oie . The mf fe range suggested by the 
electroweak data is indicated (dashed lines) for comparison. m? t ole is in GeV. 
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